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Abstract 

Using the experimental upper bound on the neutron EDM and 
experimental result on 6 ^ S7 branching ratio we have calculated CP 
asymmetry and T'^i^DM^^jj _^ sl+r)/T^'^^b sl+l'). It is shown 
that in the invariant dilepton mass q'^ region (rn-L +0.2 GeV"^) < 



PACS numbers: 12.6.Fr, ll.30.Er, 11.30.Fs 



< ni? the CP asymmetry is maximal and quite detectable. 



1 Introduction 



The experimental discovery of the inclusive and exclusive decays B Xg^ 
and B K*j by the CLEO collaboration [1,2] has triggered a lot of theo- 
retical and the experimental activity in the field of rare decays of B- mesons. 
These decays are interesting for checking the predictions of SM at one-loop 
level, for determining the CKM matrix elements, and for looking for the 
"new physics" beyond the SM. Prom the experimental point of view another 
promising decay in this direction is the semileptonic decay b — > Xgl'^l', 
because this decay is easier to measure provided that we are given a good 
electromagnetic detector and a large number of B hadrons. Theoretically this 
decay has been the subject of many works in the framework of the SM [3,4,5,6] 
and its extensions, particulary in Two Higgs Doublet Model (2HDM). 

b sl~^l~ decay is an FCNC process which appears only at the one- 
loop level of pertubation theory. The basic thing about this decay is that 
the penguin diagrams provide the two key ingredients needed for partial rate 
asymmetries. Being a loop diagram, it involves all three generations, each 
generation contributing with different elements of the CKM matrix. At the 
same time the loop effects that involve on- shell particle rescatterings provide 
the necessary absorbtive parts. 

It is well known that in 2HDM, b — > sl'^l~ decay receives significant 
contributions from the charged Higgs (H^) exchange [7]. Another interesting 
pecularity of 2HDM is the appearence of new sources of CP violation [8] in 
addition to the one in SM. An interesting version of 2HDM, so called the 
most general 2HDM, which was proposed in [9], has a new source of CP 
violation, arising from the relative phase between the vacuum expectation 
values of two Higgs scalars. 
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In this work we shall work out b — >■ sl^l~ decay. In particular we shall de- 
termine the CP asymmetry A and the ratio r = y'^hdm^Jj sl~^l~) jV^^ih —>■ 
sl^l~) as functions of the charged Higgs mass. 

In the calculation of the CP asymmetry we shall consider both the SM 
and 2HDM contributions simultaneously. In determining r and A we shall 
make use of the experimental results on BR{b — > 57) [1,2], and the neutron 
electric dilpole moment (EDM). 

Section 2 is devoted to the derivation of basic theoretical results and 
Section 3 contains the numerical analysis of them. 

2 Formalism 

In the most general 2HDM [8,9] the couplings of with tn and ba is 
characterised by the coefficients defined by 

= ■ . , e-/(^-^/) - CO0 (1) 

sinpcospsmo 

where f= t or b, (T/ = + for b and - for t, and 5f — /i2/^i where /i2 and hi are 
the diagonal elements of the matrices and P" respectively. Here P" are the 
matrices in the flavour space, and determine the Yukawa couplings (for more 
detail see [9]), and 5 is the relative phase between the vacuum expectations 
of the two Higgs scalars: 

<02> = Tf^^^'^ (2) 

The most general 2HDM reduces to the well-known 2HDM's in the current 
literature, in certain limiting cases [9] . Namely, if St = Sb = 0, then = = 
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—cotf3 (Model I) and, if = 5,5t = 0, then = —cotp, = tanjS (Model 
II). 

As mentioned above the penguin diagrams provide the necessary absorb- 
tive parts for the calculation of the CP asymmetry. In this decay the dilepton 
invariant mass ranges from Anif to ml; therefore, u and c loops give rise 
to nonzero absorbtive parts which are described, at the point /i — nib, by 

F = iAV2GFX^—AQSLiM^i^,^~ (3) 
where = VigV^l and the function Ag is given by 

Ag = w^[Q{ml/q') - Q{ml/q')] (4) 

where 

Q(x) = — (2 + 4x) VI - 4x^(1 - 4x) (5) 

9 

and Wu, having the numerical value of 0.3864, comes from the RGE movement 
of the Wilson coefficients from /i — Mw to /i = rub point. 

It is well- known that in the range (4m^, ml) one can create real low lying 
charmonium states [10,11]. In this work we shall discard that portion of total 
dilepton mass range including J/ ip and ip' poles and the region between them 
to avoid the addition of new hadronic uncertainities to the decay amplitude. 
Thus we restrict ourselves to the following kinematical regions [6] : 

Region I : Amf < q^ < (m^ — r) 
Region II : {m^, + t) < q <m^ (6) 

where r = Q.2GeV'^ is the cut- off parameter. 
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Taking into account the 2HDM contributions and absorbtive part de- 
scribed by F in (3), the amphtude for h sl^l~ can be written as 

Cio(Ai)sL7/^^i^'^7M75/~ + (7) 
— X Ci{ij)sa^i,{mbR + msL)hl^^^,l } 

The Wilson coefficients appearing in (7) are given by 

Cl^^ilj) = A,[C|^(//) + C7r^^(/i)]+a„A9 (8) 

The exphcit forms of Cf^{ij), (i=7,9,10) including leading and next-to- 
leading order QCD corrections can be found in [3,12,13,14]. The 2HDM 
contributions, Cf^^^di), in the framework of the most general 2HDM [9] 
are given by 

Cr^'^'ilA - \^t\'Kl' (9) 
where Rtb = Re[^t^t], hb = Im[itil\ and 

K - -=^^^^Biy) + yF{y) (10) 



with X = ml/M'yY-, y = m^/Mfj, s^, = 0.2315, t] = as{Mw) / as{nib) and the 
exphcit expressions for functions A, D, E, F, G can be found in [12]. 

As noted in [9], is expected to be of order of unity or less, if the Yukawa 
couphngs of the top quark is reasonable. We have shown that this happens 
to hold also for the decay process under consideration. Thus, without loosing 
generality, in what follows we set | |^ = (all the conclusions remain in 
force for the case of | |^ =1 as well). 

Using (7), the differential decay rate for h — > sl^l~ is obtained as 

j-p2HDM 9 

^-^^ = Xo{l-snA{-^+l)\CM\'+{l + 2s)i\Cf^{^l)\' + \C^M\') 
+ 12Re[CMC'/^{^)]} (11) 



where s = q'^/ml, and Aq = 

After integrating (11) over s we get 

7 = 70 + 4p/2 + 2/(6/9 + 64'^^t«) 

+ ApR^ + 2R{6R^ + 6a'ihtu + ^pC^^) (12) 

+ 1209 ^C^^Itu + o!i^ ftu + '^((J'rdhu + 0-i^Rtu) 



where 



^2HDM 



7 



Ao I At I 

rSM 
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^0 I 

I = ItbKl (13) 
R — RtbKjh 



Rtu 



I At P 
^e[AtA:] 
I AtP 
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and the other parameters in (12) are defined by the following integrals: 

,2 



P = 



,(1) _ 



/ ds{l-sf{- + l) 

J s 

R9 = J ds{l - sfRe{Cl^) 
I ds{l - sflm{Cl^) 

J ds{l - sfAg (14) 
4'^ = J ds{l-s)\l + 2s)Al 
a,9 = J ds{l - sf{l + 2s)Re{C^^)Ag 
a^g = J ds{l - sf{l + 2s)Im{C^^)Ag 

For the CP conjugate process, the analog of (12) can be obtained by the 
following replacements: 

7 = 7(7^-7; hu^-hu) (15) 

Now we introduce the parameter r that measures the relative strength of 
2HDM and SM rates 

r - ^ (16) 

ISM 

where '^sm is obtained by setting 7 = 7? = in (12). 
Next we define the CP asymmetry by 

A = t2 (17) 

7 + 7 

Substituting the expressions for 7 and ^sm into (16) we obtain a circle 
for fixed values of r: 

(7? + Rof + (7 + hf = t{r -1) + Rl + II (18) 



where the parameters Rq and Iq are given by 

lo = ^{h + a^^^Rtu) (19) 

and the quantity Tq = j^a'g^Itu is introduced for later use. 

On the other hand, insertion of (12) and (15) into (17) yields another 
circle 

{R + R!,f + + r,f ^-t + 6(1 R!^ + 7^ (20) 

where 

~ A 

K - + + l (21) 

The parameters e and t in (19) and (20) are given by 



^*^a24'^Cf^ + 2a,9; 



4p 

(1 -As) 

t = -^^47^^ (22) 
where Ag is the CP asymmetry in SM which is obtained from (17) by: 

As^A \i=i^o (23) 

Up to this point, our analysis of 6 — sl^l~ decay parallels that of 6 — 57 
in [9] except for the definition of A. We shall, however, analyze the circles in 
(18) and (20) in a different context by exploiting the relation between 7 and 
neutron EDM, and experimental results on 6 — > 57 branching ratio [1,2]. 
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First we obtain the expression for the CP asymmetry in (17) by subtract- 
ing (20) from (18) and solving for A: 

■ ^ (24) 



1-a 
where 

a = (25) 

e + 2//o + 2Rro ^ ' 

Now we turn to the determination of / with the use of the experimental 

upper bound on neutron EDM. Weinberg has proposed a CP violating 6 

dimensional gluonic operator [15] 

O, - /abcGrG^^Ge^. (26) 

which has been shown to give very large contribution to neutron EDM, d„ by 
the neutral [15] or charged [16] Higgs exchange. Weinberg, after relating the 
hadronic matrix elements of Oq to predicts the value of dn on the basis 
of a Naive Dimensional Analysis (NBA). However a detailed analysis by Bigi 
and Uraltscv [17] reports a different value for dn which equals ^ of that of 
Weinberg's. The big difference between the results of these analyses is an 
indication of the existence of hadronic unccrtainities which are mainly intro- 
duced by the matrix elements of Oq between the nucleon states. In addition 
to these theoretical unccrtainities, we have also problems with experimen- 
tal data (in that experiment yields only an upper bound on neutron EDM). 
These can be summarized as 

d'f" = CtheorXlthK{y)lQ-''''ecm (27) 

^actual ^ CexpXd^"'' (28) 

where Cthem- and c,.,xp are constants and | Cexp \ is known to be less than unity. 
Let us note that Ctheor is related to the theoretical unccrtainities and Cexp 
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to the experimental uncertainities. Experiment yields d^"^ = 1.1 10 ^^ecm 
[18]. The function K{y) in (27) is given by [16,17]: 



(29) 



The common point for the analyses in [16] and [17] is the presence of the 
function K{y) which is equal to | as y — >■ 1. 

Equating (27) to (28) and defining (3 = l.l^^a^, we obtain 

Cffi, 



eor 



I = Pf{y) (30) 

where 

/(.) = (31) 

Note that the constant f3 in (30) includes now both theoretical and exper- 
imental undeterminicies. We shall not make any assumption concerning the 
value of (3; instead we are going to fix it through the use of the experimental 
results on 6 — > S7 branching ratio. 

The 6 — > S7 decay amplitude is given by 

4:Gj? a 

M = -^—CMs{P>AmbR + msL)h{p)F^^ (32) 

where Cj{^) is defined in (8). Using the experimental result on the braching 
ratio of 6 ^ S7 decay [1,2] we get the following circle 

(Cf ^ + Rf + I^ = (Cr)2 (33) 

where is the experimental value of Crdi) 

0.22 <] |< 0.30 (34) 
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We shall determine the central values of (3, r and A which are defined in 
equations (20), (16) and (17) respectively. In doing this, we will make use 
of circles in equations (18), (20) and (33) together with equation (30). Let 
us note that (30) is obtained by the use of the experimental upper bound on 
neutron EDM [18], and (33) is constructed with the use of the experimental 
data on 6 — > S7 branching ratio [1]. 

Let us first determine (3. For this purpose we consider the circle in (33) 
in the limit of infinitely large Mh or equivalently y ^ 0. As y ^ 0, R ^ 
and through (30), / f3fo, where numerically /o = 0.2706. Then equation 
(33), which is valid for any value of Mh, yields 



(Cr)2^(Cf^^2 



p = ± r' ' ,r ' (35) 



With (35), / in (30) has now become a completely known function of M^. 
Now we solve (33) for R, yielding 

R = -C^^ + sl{CfY - P (36) 

where the choice of plus sign is necessary to satisfy asymptotic condition on 
R. 

Using (36) for R, and (30) for I we can solve equation (18) for r 

I {R + Rof + {I + Iof-Rl-Il ^37^ 

whose Mh dependence shall be discussed in the next section. 

Finally, taking r from (37), R from (36) and / from (30) we determine 
the CP asymmetry A in (24) whose dependence on Mh shall also be studied 
in the next section. 
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3 Numerical Analysis 



In the numerical analysis we shall use = lOMeV, = l.SGeV, rrib = 
4.6GeV. For the top quark mass wc rely on the CDF data [19] and for the 
W mass wc use Mw = 80.22G'e\/ [18]. 

In calculating Itu and Rtu we use the paramctrisation in [18], and in doing 
this we take the mid values of the quantities. For the phase S12 of CKM 
matrix in [18] we shall use the the mid value of cos5i^ — 0.47 ± 0.32 given in 
[20] which icludes a large uncertainity. A straightforward calculation shows 
that corresponding to the uncertainity in cos5i3, Rtu and Itu are uncertain 
by 3.87% and 23.75% respectively. Thus, the standard model asymmtery 
As in (23) is uncertain by 23.75%, and we shall use its central value in our 
calculations. This choice is justified by the closeness of Itu and Rtu calculated 
in this way to that obtained by the use of Wolfenstein paramctrisation [21]. 

Fig. 1 shows the variation of /(y) in (29) with for the lowest, central 
and the highest values of rrit permitted by the CDF data [19]. As we see 
from Fig. 1 dependence of fiy) on rrit is very weak; thus, insensitivity of 
results to the variation of y with rut is guaranteed. In what follows we shall 
use therefore the central value of CDF data rrit — 176GeV. 

For mt = 176GeV we obtain Cf^ = -0.2686. The 6 ^ 57 branching 
ratio has approximately 50% error [1] which is tranferred into a range of 
values that Cf^ may take, as described by (34). 

With the use of above-mentioned data we calculate SM CP asymmetry 
in (23) to be As = 0.0714% in Reg. 1, and As = 0.0223% in Reg. II. 

In the second column of Table 1 wc give the values of /3 as | C^^ \ moves 
from its maximum value 0.30 towards | |= 0.2686. We see that | (3 \ 

decreases gradually with decreasing | |. Moreover, it is seen that the 
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maximum value that | (3 |~ 0.5. 

Regarding the present calculations in [16] and [17] as the possible candi- 
dates for Ctheor in (27), we can make certain predictions for c^xp in (28). A 
simple calculation yields Cthear — 9.9 and Ctheor — 0.33 for Weinberg's NDA 
and Bigi-Uraltsev calculations respectively. In the case of NDA, a solution 
for Cexp exist only for | /3 |<~ 0.27 at which c?^"^*"«' turns out to be very 
close to its experimental upper bound. On the other hand, for Bigi-Uraltsev 
calculation, being a more detailed analysis, for all values of | Cj^ \ ranging 
from I C^^ I to 0.30 there exists a solution for Cex with the help of which, 
through (28), one determine the value c?"'^*""'. In the third column of Table 
1 we give the values of d'^^^^"'^ as | C^^ \ moves from its maximum value 0.30 
towards | Cf*^ |= 0.2686. We observe that for | |= 0.3 | | reaches 

its maximum value of 1.63 10~^^ which is one order of magnitude less than 
the present experimental upper bound. 

In our numerical analysis we use the range of values of Mh from AAGeV 
[18] to lOmt [15]. In Fig. 2 and Fig. 3 we show the variation of r in (37) with 
Mh in Regions / and // respectively. We observe that in both figures r is 
fairly high at low Mu and lands rapidly to a lower value after Mh ~ 500GeV. 

As wc sec from Fig. 2, dependence of r on the sign of f3 in Region / is very 
weak. Moreover, for Mh >~ ITeV, r attains the values ~ 1.056, ~ 1.0050, 
~ 1.020, and - 1.016 for p = +0.4938, -0.4938, 0.2922, and -0.2922 
respectively. 

From Fig. 3 we observe that in Region II dependence of r on the sign 
of P is large. Specificially, we see that, for large Mh, r becomes practically 
independent of Mh and attains the values ~ 1.021, ~ 0.998, ~ 1.01, and 
~ 0.9996 corresponding to /3 = +0.4938, -0.4938, 0.2922 and -0.2922 
respectively. 
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In Fig. 4 and Fig. 5 we show the variation of A in (24) with Mh in Regions 
/ and // respectively. What we observe to be common between them is the 
saturation of CP asymmetry A to a certain value after Mh ~ bOOGeV. 

Prom Pig.4 we observe that the 2HDM CP asymmetry A, practically for 
all Mh, is of the same order as the SM CP asymmetry Ag. Indeed, especially 
for large Mh, corresponding to the values of /3, /3 = +0.4938, -0.4938, 0.2922 
and -0.2922, A attains the percentage values of ~ -0.27, ~ 0.40, ~ -0.14, 
and ~ 0.28. 

In Fig. 5 we observe that asymmetry A, as compared to the previous 
figure, is completey different in that it is positive and takes higher values for 
all values of Mh- Actually, we see that for small Mh, 2HDM CP asymme- 
try is larger than the SM CP asymmetry by approximately three orders of 
magnitude. For large Mh, however, A gets values which are larger than SM 
asymmetry by two orders of magnitue. Indeed, for large Mh, corresponding 
to the values of /3, /3 = +0.4938, -0.4938, 0.2922 and -0.2922, A gets the 
following percentage values ~ 1.1, ~ 3.25, ~ 0.2, and ~ 1.5 

The last point to be noted about the Figs. 2-5 is that negative P gives 
rise to larger r and A than positive P does. 

To decern a CP asymmetry A at the a significance level with only statis- 
tical errors, the number of B hadrons Nb needed to demonstrate the asym- 
metry is given by[22] 

Now denoting the number of B hadrons to observe Ag, A in I and A in II 
by Nh, Nh and respectively, we get, using the values of r and A we have 
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obtained already, the following ratios 



10-* (39) 



which clearly prove that Region II is more suitable for experimental inves- 
tigations on A. 

In conclusion we have determined the 2HDM CP asymmetry A, ratio of 
2HDM decay rate to SM decay rate r and actual value of neutron EDM. In 
doing these we have utilized the experimental results on b ^ 57 branching 
ratio, and on the upper bound of neutron EDM. Both r and A relax to 
constant values after ~ 500GeV. This saturation property of quantities 
shows that if charged Higgs mass happens to be large (~ ITeV) then the 
most general 2HDM merely shifts the SM values of r and A to some other 
value which may be important for establishing 2HDM. Boldly speaking, in 
the high dilepton mass region (Region 77) r is closer to unity and asymmetry 
is very large as compared to those in low dilepton mass region (Region 7). 
Thus on the basis of the order of magnitude analysis carried out for Nb, we 
conclude that the high dilepton mass region is important and appropriate 
for experimental check of the quantities under concern. Region 77 [6] is 
accessible to the B experiments which will be carried out with hadron beams 
in CDF, HERA and LHC. 
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Figure Captions 

Figure 1: The Mh dependence of f{y) for rrit — 194:GeV (with circles ), 
rrit — 176GeV (bare sohd curve) and rrit — 158GeV (with squares). 

Figure 2: The Mh dependence of r in Region I. Here labes 1, 2, 3 and 4 
correspond to /3 = 0.4938, -0.4938, 0.2922 and -0.2922 respectively. 

Figure 3: The same as in Fig. 2 but for Region //. 

Figure 4: The dependence of A in Region /. Labels have the same 
meaning as in Fig.l. Here the unlabled solid line shows the SM asym- 
metry. 

Figure 5: The same as in Fig. 4 but for Region II. 
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